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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  MEMORANDUM  1334 


TEE  EFFECT  OF  HIGH  VISCOSITY  ON  THE  FLOW  AROUND 
A  CYLINDER  AND  AROUND  A  SPHERE*^ 
By  F.  Homann 


For  the  determination  of  the  flow  velocity  one  is  accustomed  to 
measure  the  impact  pressure,  i.e.,  the  pressure  intensity  in  front  of 

an  obstacle.   In  incompressible  fluids  the  impact  pressure  is  /v  /2g 

\y }      ^S,/'^  }    is  specific  weight  and  v,   m/sec,  is  velocity)  if  the 
influence  of  viscosity  can  be  neglected.   Such  an  influence  is  appreci- 
able, however,  when  the  Reynolds  number  corresponding  to  impact  tube 
radius  is  under  about  100,  and  must  consequently  be  considered,  if  the 
velocity  determination  is  not  to  be  faulty.   The  first  investigations 

of  this  influence  are  included  in  the  work  of  Miss  M.  Barker  .   In  the 
following  pages,  experiments  will  be  reported  which  determine  the  inten- 
sity of  impact  pressure  on  cylinders  and  spheres;  furthermore  a  theory 
of  the  phenomenon  will  be  developed  which  is  in  good  agreement  with  the 
measurements . 

The  research  apparatus  consists  of  an  oil  circulation  in  which  the 
velocity  of  the  oil  can  be  varied  from  0.5  centimeter  per  second  to 
30  centimeters  per  second  with  the  help  of  a  vsine-type  pump  lying 
entirely  in  the  oil.   A  Russian  bearing  oil  and  a  mixture  of  this  with 
fuel  oil  is  used  for  the  measurements.   Figure  1  illustrates  the  test 
setup.   In  this  is  indicated:   P,  the  pump;   a,  turning  vanes; 
G,  straightener;  and  V,  the  actual  test  section  which  possesses  a 
breadth  of  0.148  meter,  a  depth  of  0.15  meter  calculated  from  the  oil 
surface,  and  a  length  of  O.Jk   meter.   It  was  provided  with  wall  ports  A 
in  three  different  places.   E  is  an  entrance  section  for  the  pump; 
D,  a  diffuser;  the  immersion  heater  T  and  the  cooling  coil  K  provide 

"Der  Einfluss  grosser  Zahigkeit  bei  der  Stromung  um  den  Zylinder 
und  um  die  Kugel",  ZAMM,  vol.  l6,  no.  3,  June  I936,  p.  153-164. 

-^The  suggestion  of  the  present  work,  which  was  prepared  in  the 
Kaiser  Wilhelm  Institute  in  Gottingen,  I  obtained  from  Herr  Professor 
Dr.  Prandtl,  to  whom  in  this  place  I  express  my  most  heartfelt  thanks 
for  the  energetic  furthering  of  the  work  and  the  valuable  suggestions 
given  me  for  its  completion.  Another  work  in  the  same  field  is  published 
in  "Forshung  a.d.  Gebiet  des  Ingenieurwesens",  I936,  vol.  7,  no.  1. 

T^.  Barker  -  Proc.  of  Royal  Society,  1922,  Vol.  A. -101,  p.  4-35. 
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temperature  regulation.   For  impact  pressure  measurement  a  cylinder 
which  was  provided  with  a  port  was  built  rigidly  into  the  test  section. 
The  dieuneters  of  the  cylinders  used  were  1  centimeter,  1.377  centimeters, 
1-953  centimeters,  and  2 'centimeters.   The  holes  had  a  diameter  of 
0.1  centimeter  and  0.2  centimeter.   Two  corrections  -  one  because  of 
wall  effect,  the  other  because  of  finite  size  of  hole  (which  originated 
with  A.  Thorn  and  first  had  to  be  checked  for  the  measuring  range  under 
consideration  here)  -  were  applied  to  the  measurements,  which  are  illus- 
trated in  figure  7.   The  solid  curve  represents  the  theory. 

More  precise  information  on  the  test  setup  and  the  measurement 
technique  are  found  in  the  work  cited  in  footnote  1. 

In  the  case  of  the  measurement  of  static  pressure  on  a  sphere,  a 
sphere  provided  with  a  hole  was  affixed  to  a  pitot  tube,  the  sphere 
having  one  time  a  diameter  of  O.B  centimeter,  the  other  time  1.6  centi- 
meters.  The  execution  of  the  measurements  was  in  the  same  manner  as  in 
the  case  of  measurement  of  the  Barker  effect.   The  result  is  shown  in 
figure  8.   The  solid  curve  corresponds  to  the  theory. 

In  order  to  arrive  at  a  clearer  picture  of  the  viscous  flow  around 

a  cylinder  or  a  sphere,  the  case  of  viscous  flow  against  a  plate  was 

next  calculated.   The  differential  equation  appearing  in  the  two- 

■:> 
dimensional  case  has  already  been  solved  by  Hiemenz-^  and  will  be  sketched 

once  more  for  the  sake  of  a  better  understanding  of  the  final  form.   The 

solution  will  than  be  used  in  the  flow  around  the  cylinder. 

After  this,  the  three-dimensional  flow  against  a  plate  in  a  fluid 
jet  will  be  treated,  to  be  used  on  the  sphere. 


VISCOUS  FLOW  IN  THE  VICINITY  OF  A  STAGNATION  POINT 
(TWO-DIMENSIONAL  CASE) 


Potential-flow  theory  gives  for  the  velocity  components  in  the 
neighborhood  of  a  stagnation  point,  for  the  case  of  flow  perpendicular 
to  a  plane  wall  (fig.  2): 


U  =  -ax     V  =  ay 


-%iemenz  -  Dissertation,   Dingl.    Polyt.    Journal,   V.326   (I9II), 
No.    21-26. 
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The  pressure  is  found  from  the  Bernoulli  equation  to  be 

,2 


P^  -  P   2  -   ■  ■  -    2 


a(u2  ,  v^)  =  P|l(x2  .  y2) 

With  consideration  of  viscosity,  Hiemenz^  makes  the  formulation: 

u  =  _f(x)     V  =  yf'(x)  (1) 

^o  -  P  =  ^(^x)  ^  ^)  (  =  ) 

The  continuity  equation  is  fulfilled;  the  boundary  conditions  read: 

for  X  =  0  (that  is,  at  the  wall):   u=v=0,   f=f'=0 
for  X  =  co;   V  =  V,   f  =  a 

In  equation  (2)  if  p^   signifies  the  pressure  at  the  stagnation  point, 

then  F/  v  =  0.   From  the  equations  of  motion 
(o) 

sx       ay      p  ax      \;d^^    ^/l 
^-       ^y     p^y      W    ^yV 

one  obtains  as  determining  equations  for  f,    x   gind  F/  x  : 

(x)        (x) 

a2 

ff«  =  f_  F'  -  Vf'»  (3) 

f»2  _  ffit  =  a2  +  vf"  (k) 

with  the  boundary  conditions  above. 
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Equation  (k)   has  already  been  integrated  by  Heimenz.   In  order  to 
make  the  coefficients  ej^ual  to  unity,  he  sets 

f(x)  =  A^(0     I  =  ax 

From  comparison  of  the  coefficients,  it  follows  that 

A  =  \A^     a  =^ 

With  this,  equation  (4)  becomes 

0t2  _  ^,,^  ^  -L  ^  0,,,  (5) 

The  new  boundary  conditions  read 

for   |=O:^=0'=O 
for   ^  =  00:   0»  =  1 

The  behavior  of  0  and  the  first  two  derivatives  is  shown  in  figure  4. 

We  need  the  pressure  difference  between  the  stagnation  point  and 
the  pressure  for  x  =  ocw   For  x  =  00: 

^'  =  1     f  =  a     5?  =  I  -  0.647 

From  integration  of  equation  (3)^   F  is  determined  to  be 

^  F  =  Vf '  +  I  f 2  (6) 

If  one  now  forms   (p  -  p^  minus   (Pq'  -  P')^  ^.s  given  by  the  Bernoulli 
equation,  one  obtains: 

(P„  -  P)  -  (P,-  -  P')  =  if  (r(,,  .  y^)  -  ^-{fj  *  y%„.^) 
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For  the  stagnation  streamline,  for  which  y  =  0: 

(p„-p)-(p„'-P')=ef^F(^,  |i-.2 

If  one  puts  in  for  F  the  previously  obtained  value,  there  results 

(P^  -  P)  -  (Pq'  -  P')  =  |(2Vfeo'  +  ^I    -    fco^)  =  PVfeo' 

For  X  =  00,   f •  =  a;  therefore,  one  obtains  as  a  final  formula: 

(Pq  -  p)  -  (Pq'  -  P')  =  Pva  (6a) 

VISCOUS  FLOW  AT  A  STAGNATION  POINT 
(ROTATIONALLY  STOMETRIC  CASE) 


For  the  solution  of  the  differential  equation  arising,  all  the 
expressions,  such  as  the  equations  of  motion,  the  velocity  components, 
etc.,  were  reduced  to  cylindrical  coordinates. 

If  z,   r  and  |3  are  the  coordinates  (fig.  3)^  then  corresponding 
to  the  two-dimensional  case,  there  will  apply: 


V.  =  -f/„^     v^  =  J  f(^)      .  (7) 


P„-P 


The  continuity  equation  is  again  fulfilled;   v„  =  0,  since  we  are  dealing 

with  a  rotationally  symmetric  process.   These  expressions  stem  from  the 
frictionless  problem  of  a  fluid  jet  against  a  plate,  where 

Vj.  =  ar     72  =  -2az 

and 

p   _  p  =  P|!(i,,2  ,  ,2) 
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The  quantity  2az  in  the  frictionless  case  is  replaced  by  f/  x   in 
the  viscous  case.   In  the  case  at  hand  the  equations  of  motion  read: 


dVj,       dvj, 
or       dz 


Bv„      Sv„ 
^  dr     ^  dz 


=  .1  ^  + 
P  Sr 


^a^v 


-^    dVv,    V 


r   X  ^'  r 
^  dr 


r2    az^y 


,1  ^ 
P  Sz 


^S^v 


av„     a^v. 


+  V 


^  +  11:^  + 


^az 


2    r  Sj 


az^ 


>    (9) 


Substituting  equations  (7)  and  (8)  in  equation  (9)  gives; 


i  ft2  _  f-ftf  =  2a2  +  vf*'* 


f  f  •  =  §:_  F»  _  vf ' ' ' 


(10) 
(11) 


The  boundary  conditions  read 


for  z  =  0:   f  =  f '  =  0 
for  z  =  00:   f  =  2a 


If  one  finds  f  from  equation  (lO),  one  can  therewith  determine  F 
from  equation  (ll).   One  next  substitutes  the  transformation 


f(z)  =  A0(O     %   =  az 


(12) 


into  equation  (10),  in  order  to  make  the  two  coefficients  equal  to  unity. 
This  yields: 

i  a^A^0'^  -  a^A^00"  =  2a2  +  va^A^'" 


From  equating  the  coefficients; 


12  2     2     3 
-  a  A  =  2a  =  Va-^A 
2 
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A  =  2^     ^  =^  (13) 

From  equation  (lO)  vrith  equations  (12)  and  equations  (I3)  there  results 
the  final  differential  equation: 

(3»'«  +  25;>0"  -  52'^  +  1  =  0  ilk) 

with  the  boundary  conditions: 

I 

for   1=0:   0  =  0'  =  0 
for  I  =  00:  0«  =  1 

The  differential  equation  (l4),  just  as  Hiemenz',  is  no  longer 
elementarily  integrable.   Its  solution  was  obtained,  accordingly, 
through  a  power  series  development  from  zero: 

^  =  a^  +  a^^l  +  agl   +  .  .  .  +  a^^^  (15) 


By  the  method  of  undetermined  coefficients,   a^^   can  be  determined. 

Since,  however,  one  boundary  condition  lies  at  infinity,  one  coefficient 
remains  undetermined;  and  in  fact  it  turns  out  to  be  ao.   From  the 

recursion  formulas 

2      3 
P  =  a.^  +   a-j^l  +  a  I  +  a^l"^  +  .  .  . 

0'  =  a  +  2a2|  +  3a^|  +  .  .  . 

0"  =  2a2  +  2  X  3a3l  +  3  X  4a.  |^  +  .  .  . 


result  as  coefficients: 
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a^   =  a^  =  0 

a2  =  for  the  present  undetermined 

a^   =  -O.I6666T 

a^  =  0 

a^  =  0 

a^  =  0.555556  X  10" 2  ag 

a^  =  -0.396825  X  10" 3 
ag   =  0 

a^  =  -0.41+0917  X  10" 3  a.^ 
^10  =  0.793651  X  10"^  ag 
^11  =  -0.360750  X  10-5 
^12  =  0.371+111  X  10"^  B.^ 
a-j^2  =  -0.114597  X  10"^  a.^ 
a-^^  =  0.115735  X  10-5  ag 
/  ai5  =  -0.301lti32  X  10"5  B.^  -  0.385784  X  10"''' 

a^g  =  O.I34S96  X  10~5  agS 
a^      =  -0.211005  X  10"^  ag^ 

^18  "  0.224141  X  10"^  a25  +  0.157758  X  10"'''  ag 
a-j^^  =  -0.135546  X  10"^  ag^  -   0.415153  x  lO'^ 
ago  =  O.3I6633  X  10"'''  ag3 

agi  =  -0.152798  X  lO-T  ag^  -   0.295658  X  10"^  ag^ 
agg  =  0.119505  X  10"'''  ag5  +  0.199390  X  10"9  ag 
ag3   =  -0.371665  X  10"^  ag^  -   0.433457  X  10"^! 
agi^   =  0.956242  X  10-9  a^T  +  0.55^360  X   10-9  ag3 
a25  =  -0.943031  X  10-9  ag6  -   0.462914  x  10"^^  ag2 
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In  order  now  to  be  able  to  determine  ag^  a  second  series  develop- 
ment from  infinity  was  set  up,  which  was  adjusted  to  the  boundary  condi- 
tion for  0  at  infinity.   To  this  end  one  sets 

■  0  =  0o  +  S^i  (16) 

in  which  0-,  corresponds  to  a  small  quantity,  which  one  can  neglect  in 
the  following  expressions  when  it  appears  squared.  ^^  is  the  solution 
for  ^  =  00. 

0.  =  0^'  +  0-^'  =  1  +  0i' 

since  for  |  =  oo; 

0^t  =  0t  =  1     0tt  ^  ^^tt     0tfi  =  ^^...         (i6a) 

The  boundary  condition  reads 

for   I  =  oo:   0  «  =  0 
Furthermore,   0  =  |. 

The  integration  constant  is  omitted,  since  in  the  following  calculation 
it  comes  in  again  automatically. 

If  one  substitutes  the  above  values  into  equation  (lO),  one  obtains 

i^l'"  +  2(^0^1"  ^  ^A")  -  (^o'^  ""  ^o'^l'  +  ^1'^)  +  1  =  0  (IT) 
Or  if  one  neglects  the  squared  terms  in  0  : 

0^'«'  +  2|0^"  -  20^'  =  0  (18) 


10  NACA  TM  1334 

To  solve  this  differential  equation  one  sets 

$  =  0i' 

With  this,  equation  (I8)  gives 

$f t  +  2|$»  _  2$  =  0  (19) 

A  special,  not  identically  vanishing  solution  of  equation  (I9)  is: 

If  Oo  is  an  additional  solution,  then 


21  dl 


'^1^2'  "  ^l'^2  ~  ^ 


_p2 


*2    |2|^      ^ 


This  equation  is  directly  solvable.   Its  solution  is 

^00   T]^ 


The  general  solution  of  equation  (I9)  is  then: 

00  T 


«D  =  c^o-L  +  CgOg  =  c^i  +  ^2^1     ^  ^"^^  'in 
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11 


Since   for     ^  =  00,      J^   '    =   $  =  0,    then     C-j_   =  0.      Therefore 


0/    =  0   =  C2I  /       -^  e-n^  dTi   =  Cgl 

-'-  J  00       TJ 


ie-^"-2/'e-1^d, 


and   since 


<^-^   =  *  d? 


03_  =  C^   +  Cg  j 


2  pn      ,,2 

■e-n    -  2^  e 

<-'  00 


-^     du 


dn 


I 


Ml  2 


C3   -  Cp    /        e-T     dT]   -   2C2    /        Ti  dTi    /       e-^     du  (20) 


00  ^00 


The  double  integral  becomes,  according  to  Blasius  \, 


2C. 


T\  d^     /   e"^  du  =  2Cr 


1  «=2 
2 


''L  ^'''"^-il  --^' ^^  *  i  ^--^^ 


With  this,  equation  (20)  becomes: 


01  -  C 


I  2 


Sr^-M-^  -ft^^^)/  -^ 


dTl 


(21) 


Now: 


o|      2  r^      2  r~      2  r^      2 

2   /        e"^     dT)  =   2   /        e-^     dT)  -   2   /        e"^     dr]  =   2   /        e'T     dr^ 
^00  ^o  ^o  ^o 


^ 


\AtJ 


r  I        2 
e   ■•     dT)  =  -f=  /        e'T     dTl   -   1 


>(22) 


^. 


o 


Blasius  -  Dissertation,  Z.F.M.u.  Phys.  1908^  P-  !• 
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If  one  substitutes  equation  (22)  in  equation  (21),  one  can  calculate  (f> 
pointwise,  since 


2  r^  -.2 


\/s, 


e-11  dT] 


o 


is  tabulated.   Therefore 


^.-i-^^-(i-^)f{i/;-^^-: 


t=^[-i^x''="'''4 


Herewith  0,   0',  and  0*'   are  determined  for  the  development  at 
infinity,  as  a  comparison  with  equations  (l6)  and  equation  (l6a)  shows. 

In  both  developments   ao,   C2,  and  Cn   appear  as  unknowns.   If 

one  now  combines  both  solutions  at  the  point  \=  \q,    and  determines 

that  the  value  of  the  function  and  the  first  two  derivatives  of  the 
series  development  at  zero  are  equal  to  the  corresponding  values  that 
one  obtains  from  the  development  at  infinity,  then  three  determining 
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Thereby  ao  is  determined  accurately  to  at  least  five  places.   For  the 
coefficients  of  the  power  series  this  yields: 


a„  =  0.658619 


^2 

^3  = 

%  = 

"9  = 

^10  = 

^11  = 

^12  = 

^13  = 

hk  = 


-0.166667 
0.365900  X  10"^ 
-0.396825  X  10" 3 
-0.191261  X  10" 3 
0.52271^  X  10"^ 
-0.360750  X  10" 5 
0.106882  X  10"^ 
-0.497098  X  10"5 
0.762253  X  10"^ 


^16 
^17 
^18 
^19 
^20 
^21 
^22 

"23 

^2k 

"25 


0.385391  X  10-6 
-0.958673  X  10"''' 
0.381678  X  10"''' 

-0.259200  X  10"'^ 
0.904605  X  10"^ 
-0.252966  X  10"° 
0.161233  X  10"^ 
-0.703675  X  10"9 

0.209783  X  10"^ 
-0.970520  X  10"^° 


a^  =  -0.605859  X  10 


-6 


a^  -  aj^  -  a^  -  Bg  -  0 


The  values  for  ^,      0',  and  0''   are  shown  more  accurately  how- 
ever, in  Table  I.   In  this  case,  0''   is  calculated  accurately  to  two 
decimal  places,   0'   to  two,  and  0  to  three.   With  this  the  differ- 
ential equation  (l4)  is  solved. 

From  integration  of  equation  (ll),  one  obtains 


3^  F  =  vf '  +  i  f 2  =  2av(0'  +  0^) 


(23) 
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As  in  the  plane  case,  one  uses  again  the  pressure  difference 
between  the  stagnation  pressure  and  the  pressure  for  z  =  00.   For  |  =  0, 
equation  (23)  is  equal  to  zero;  for  |  =  00 

0'  =  1     f  =  2a     0=1-  0.557611 

If  one  now  forms  again  /p^  -  p\   minus   (p^'  -  p'),  as  given  by 
the  Bernoulli  equation,  one  obtains: 

(Po  -  p)  -  (Po'  -  p*)  =  p^^"'  -^  i  ^~'  - 1  ^-')  =  p^^»' 

As  a  final  formula  one  obtains 

(Pq  -  p)  -  (Pq*  -  P')  =  2pva  (2I+) 

STAGNATION  PRESSUEE  ON  A  CYLINDER 


For  the  stagnation  streamline  the  Navier-Stokes  differential  equa- 
tion gives 

u^  +  i^=v^+v^  (25) 

Figure  12  shows  the  variation  of  u  on  this  streamline.   The  different 
behavior  of  u  at  the  stagnation  point  from  potential  flow  is  explained 
by  the  influence  of  viscosity.   If  one  integrates  between  the  boundary  R 
and  00^  one  gets 


"k" 


"^^^pH-M-ri^-rg" 


Figure  I3  shows  that 


V  ^ 


R 

=  0 
00 
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likewise  Uj^  =  0,  and  we  want  to  identify  pp  with  p^   of  the  pre- 
ceding calculation.   One  obtains,  therefore 

dU  2      r^  5.2 
(Po  -  P)  -  ^  =  P^  /   H  ''^  ^^^^ 


00   Sy 


■  R  2 

2-ii  dx 

°  Sy2 


is  calculated  approximately  in  that  for  u  the  value  corresponding  to 
the  potential  flow  is  put  in.   The  contribution  of  the  boundary  layer 
to  the  integral  is,  in  the  case  of  not  too  small  Reynolds  number,  small 
in  comparison.   As  potential  function  $   of  the  flow  around  the  cylinder, 
one  obtains 


=  4  ^  f. 


and  with  it: 


^  =  _u  Z'-  25!  +  ^^  +  8y?R!  _  ^Qx^ygRg' 


For  y  =  0,  therefore,  along  the  stagnation  streamline: 


A  2      h 
ay  X* 


-r&-=^ 


Herewith  equation  (26)  gives 

(Po-P)=^-^  (27a) 


p\]J        2pvU, 
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or 


^^  =  1  .  ^  (27) 

pU„2/2       ^- 


If  one  substitutes  7  =  pg,  then  formula  (2?)  reads 

7UoV2g     R^ 

where  g  =  9-8l  meters  per  second  . 

In  order  to  be  able  to  accomplish  a  comparison  of  test  results  with 
theory^  the  "displacement  thickness"  (see  Tolmien:  Hdb.  d  Experimental- 
physik,  V.  k,    1st  part,  p.  262,  "Grenzschichttheorie  (Boundary  Layer 
Theory)")  on  the  cylinder  must  yet  be  considered  in  the  calculation. 
Solution  of  the  differential  equation  (5)  yields  (fig.  5): 

I*  =  ax*  =  a6*  =  0.6ij-7 
where  5*  is  the  displacement  thickness.   Therefore 

0.647  =  6*/|  (28) 

If  one  compares  the  flow  in  the  region  nearest  the  stagnation  point  for 
the  cylinder  and  for  the  flow  against  a  plate,  one  obtains  from  equa- 
tions (27a)  and  (6a) 

2pvUo       2Uo 
pva  =  -^,      a  =  -^ 

If  one  substitutes  this  value  in  equation  (28),  one  obtains  for  the 
displacement  thickness 

0.61+7  =  \/^  S*  =  ^\^Ri 
VRV       R 

The  dependence  of  6*/R  on  Re  is  shown  in  figure  6. 
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Now  since  in  the  test  results  Re  is  formed  from  the  cylinder 
radius  R,  the  actual  effective  radius  is  therefore   (R  +  5*),  and  equa- 
tion (27b)  is  altered  to: 

^°  "  ^  =  1  +     ^^ 


„  2/_       Uo(R  +  5*) 
7U0  /2g        °^  ' 


With  this  one  obtains  as  a  final  rule  for  the  stagnation  pressure  on  the 
cylinder 

"°  ■  '-  =  1  .  ^ (29) 


7Uo^/2g       Re  +  0.il57^ 


In  figure  7  "the  solid  curve  again  gives  the  theory,  which  agrees  very 
well  with  the  practice. 


STAGNATION  PRESSURE  ON  A  SPHERE 


Corresponding  to  a  cylinder,  for  the  stagnation  streamline  of  a 
sphere 


au     1  ap       a^u      /a^^     s2u^ 

u  +  —  -^^  =  V  +  v( +  

ax   P  dx     ^^  Uy2   ^^2 


If  one  integrates  again  over  x  from  00  to  R: 


(Po-P)^^-v/;(0.^)a,  ,30, 


Since 


.R  .2 

o  u 


dx  =  0,   %  =  0 


^  00  Sx*^ 
The  integral  on  the  right  side  of  equation  (30)  one  again  solves  by 
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substituting  for  u  the  value  for  the  potential  flow.   The  potential 
function  of  this  flow  is 


* '- '""('  *  S3) 


For  potential  flow  it  is  further  true  that 

2  2? 

S  u        d  u        9  u        - 

+  — 7=;  +  — ^  =  0 


Therefore 


ox         dy         oz 


(^  u       S  uK     _    r      B  u    ,  5u 


R 


If  one   substitutes   in  the  last   formula  the  value   for     bu/bx     given  by     0, 
one   obtains 

m  ■  »■  ■  ^ 

Herewith  equation  (30)  becomes 


2 

)_ 
'o   "     2       R 


p„  -  p  .  £^  .  3PVU,  ,^^^j 


or 


^°  "  '  =  1  +.  ;^  (31) 


pUoV^ 


Re 


or  with  7  =  pg 


''-'    =1^4  (31.) 


7u//2g 


Re 
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If  here  one  also  puts  the  displacement  thickness  Into  the  calculation, 
one  gets  (since  in  the  rotationally  symmetric  case  I*  =  0.5!)76): 


0.5576  =  &*Ja/^ 


Comparison  of  equation  {2k)   and  equation  (jla)  yields 

.,  -  3Uo 


2R 


0.5576  =  \/^  B* 
V  2VR 


5^  ^  0.^55 
R 


\fRe 


(32) 


It  appears  that  the  displacement  thickness  I'or  a  sphere  and  a  cylinder 

are  equal  within  ^  percent,  although  the  displacement  thickness  in  the 

case  of  plane  flow  against  a  plate  is  different  from  the  corresponding 
three-dimensional  flow. 

If  one  considers  the  displacement  thickness  in  equation  (31b),  one 
obtains  as  a  final  stagnation  pressure  formula  for  a  sphere 


P     -   P 
^o        ^ 


=  1  + 


Re  +  0.it55v/Se 


(33) 


The  solid  curve  in  figure  8  corresponds  to  the  theory;  the  agreement 
with  test  results  is  again  satisfactory. 

From  the  final  stagnation  pressure  formula  the  dependence  of  the 
numerical  factor  c   on  Re  can  be  determined,  if  one  sets 

Pq  -  P  ^  1  ^  _!. 
7Uo72e  '   ^  ^^ 

For  the  sphere  there  results 


e  = 


6  Re 


Re  +  0.i+55v^ 


(3^) 
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From  Stokes'  calculation  one  obtains  for  small  Re:   e  =  3.   In 
figure  9  is  drawn  log  Re  as  abscissa,   g   as  ordinate.   In  the  region 
from  about  Re  =  0.1  to  Re  =  1  the  course  of  c   is  essentially  dif- 
ferent, since  Stokes'  law  describes  an  approximation  for  very  small 
Reynolds  number  and  the  above  law  is  an  approximation  for  large  Reynolds 
number. 

For  the  cylinder  one  obtains  in  the  same  fashion 

e  =  ^^V-^  •      (35) 

Re  +  O.h^lVRe 

According  to  Lamb-^,  for  small  Re,  for  which  the  validity  of  the  formula 
extends  to  about  Re  =  0.5: 

^^_^e_^ Re  (36) 

1.309  -  ZnRe 


In  figure  10  is  again  shown  the  dependence  of  g   on  log  Re.  Within 
the  accuracy  of  measurement  the  test  results  here  also  confirm  equa- 
tion (29). 

With  the  help  of  the  flow  against  a  plate  it  is  now  also  possible 
to  establish  approximately  the  course  of  u,   du/Sx,  and  from  this  p, 
on  the  stagnation  streamline.   A  single  curve  was  assumed  in  which, 
inside  the  displacement  thickness  &*,  the  magnitudes  as  given  by  the 
flow  against  a  plate  were  used.   From  the  displacement  thickness  on, 
which  had  a  value  of  O.Oil55  in  the  foregoing  case  for  Re  =  100,  the 
potential  flow  was  calculated.   To  explain  the  transition  from  viscous 
to  potential  flow,  I  would  like  to  go  through  the  calculation  of  u  as 
an  example.   The  solution  of  the  viscous  problem  u,  /u   has  as  asymptote 

the  tangent  to  the  curve  U2/uq,  which  was  determined  from  potential 

theory,  at  the  point  6*  =  O.Oij-55  centimeter.   In  figure  11  this  tangent 
is  labelled  t.   The  difference  k  between  the  asymptote  t  and  u-]_ 

at  the  point  x^  gives  the  deviation  of  viscous  flow  from  potential 

flow  at  this  point.   Therefore  to  the  value  u-.   at  the  point  x^  was 

added  the  proper  k.   With  the  help  of  this  procedure  one  obtains  point- 
wise  the  transition  from  Un   to  uo. 


-^amb  -  "Hydrodynamics"  (2nd  Edition  1931;   German  Edition  by 
Helly,  p.  696,  par.  3^3). 
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Su/3x  was  determined  correspondingly;  the  pressure  p  was  found 
from  the  equations  of  motion  to  be,  in  the  case  of  the  sphere: 

^  "^  =1-''A-^^^-  ^—^  (37) 


7UoV2g       ^00   200  ax    Uo(R  +  x)^ 


Instead  of  6  in  the  last  term  of  the  preceding  equation  (37)^  in 

the  case  of  the  cylinder  one  gets  the  factor  k.      Figures  12,  13,  and  ik 

are  the  results;  by  way  of  comparison  the  corresponding  curves  for  the 

cylinder  and  the  sphere  are  shown  on  one  sheet.   The  curves  are  true, 

as  already  said,  for  Re  =  100,  in  which  R  =  0.01  meter;  Uq  =  1  meter; 

2 
V  =  0.0001  kilogram  X  second  per  meter   was  assumed. 

Naturally  the  last  curves  give  only  an  approximation,  which  can  be 
made  essentially  better  through  a  second  approximation;  yet  this  task 
in  the  framework  of  the  foregoing  work  would  lead  too  far. 


SUMMARY 


In  the  foregoing  work  the  stagnation  pressure  increase  on  cylinders 
and  spheres  brought  about  through  the  influence  of  large  viscosity,  was 
reported  on. 

For  the  three-dimensional  problem,  hence  the  flow  around  a  sphere, 
a  differential  equation  was  set  up  which  corresponded  to  that  of  Heimenz, 
who  had  already  solved  the  two-dimensional  case.   The  solution  was 
ascertained  likewise  through  an  approximate  method.   The  solutions  for 
the  two-  and  for  the  three-dimensional  case  were  used  for  the  flow 
around  the  cylinder  and  sphere  respectively;  the  formulas  so  obtained 
for  the  stagnation  pressure  increase  stood  in  good  agreement  with  the 
reported  test  results.   Finally,  a  procedure  to  determine  the  velocity 
and  pressure  variation,  as  well  as  the  variation  of  bu/bx     on  the 
stagnation  streamline  was  shown  and  used  on  the  practical  case  of 
Re  =  100. 


Translated  by  D.  C.  Ipsen 
University  of  California 
Berkeley,  California 
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TABLE  I 


1 

0 

0. 

0" 

1 

0 

0. 

r 

0 

0 

0 

1.3172 

1.1+ 

0.851+6 

0. 91+76 

0  . 1697 

.1 

.006  U 

.1267 

1.2172 

1.5 

.9502 

.9635 

.1301 

,2 

.0250 

.2434 

1.1173 

1.6 

1.01+72 

.9762 

.0895 

.3 

.05J+8 

.3502 

1.0181 

1.7 

1.1^^53 

.9863 

.0622 

.k 

.097^ 

.l^U71 

.9200 

1.8 

I.2I+2I+ 

.9905 

.01+18 

.5 

.ii^39 

.53i^3 

.8235 

1.9 

1.31^36 

.9935 

.0276 

.6 

.2012 

.6129 

.7298 

2.0 

1.1+1+30 

.9962 

.0180 

.7 

.2659 

.680I1 

.6U00 

2.1 

1.51+13 

.9979 

.0115 

.8 

.3370 

.7i^00 

.55i+8 

2.2 

I.6I+09 

.9986 

.0073 

.9 

.^137 

.781+7 

.h7k2 

2.3 

1.71^16 

.9991 

.001+2 

1.0 

.U951 

.8352 

.U015 

2.1+ 

1.81+17 

.9995 

.0027 

1.1 

.5805 

.8712 

.3351 

2.5 

1.9I+20 

.9997 

.0016 

1.2 

.6653 

.9025 

.2760 

2.6 

2.0I+23 

.9999 

.0009 

1.3 

.7608 

.92I17 

.221+1 

2k 


NACA  TM  1334 


Figure  1.-    Test  tunnel. 


Figure  2.-    Streamline  picture  for  flow  against  a  plate  (two-dimensional). 


Figure  3. 
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Figure  4.-    0,  0',  0".    The  curves  drawn  out  illustrate  the  two-dimensional 
solution,  those  not  drawn  out  the  three-dimensional. 
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Figure  6.  -    Momentum  thickness  on  a  cylinder    —  =  MiJ.  and 
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Figure  7.-    Stagnation  pressiire  on  cylinder. 
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Figure  8.-    Stagnation  pressure  on  sphere. 
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o       Test     Point 
Curve   i:  £  =  3  (stokes   solution) 


Curve  n:€ 


gRg 


-3 


Re*o.45  5\jRe 


8 

o 

0           " 

^  T 

■S^^o^ 

i-^ 

X 

y 

^ 

-2 


0  I  2 

^Log    Re 


Figure  9.-     €  for  sphere. 
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Figure  10.-    e   for  cylinder. 
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Figure  11.-    Illustration  of  interpolation  for  transition  from  viscous  to 

potential  flow. 
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Figure  12.-    Velocity  variation  on  stagnation  streamline;  cylinder:    Curve  I; 

sphere:    Curve  II.    Re  =100. 
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Figure  13.-    Variation  of    —    on  stagnation  streamline;  cylinder:    Curve  I; 

sphere:    Curve  U.    Re  -  100. 
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Figure  14.-    Pressure  variation  on  stagnation  streamline;  cylinder:    Curve  I; 

sphere:    Curve  H.    Re  -  100. 
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